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Abstract. Based on a detailed study of the temperature structure of the intracluster medium in the halo of M 87, 
abundance profiles of 7 elements, O, Mg, Si, S, Ar, Ca, and Fe are derived. In addition, abundance ratios are 
derived from the ratios of line strengths, whose temperature dependences are small within the temperature range 
of the ICM of M 87. The abundances of Si, S, Ar, Ca and Fe show strong decreasing gradients outside 2' and 
become nearly constant within the radius at ~ 1.5 solar. The Fe/Si ratio is determined to be 0.9 solar with no 
radial gradient. In contrast, the O abundance is less than a half of the Si abundance at the center and has a 
flatter gradient. The Mg abundance is '-^l solar within 2', which is close to stellar abundance within the same 
radius. The O/Si/Fe pattern of M 87 is located at the simple extension of that of Galactic stars. The observed 
Mg/0 ratio is about 1.25 solar, which is also the same ratio as for Galactic stars. The O/Si/Fe ratio indicates 
that the SN la contribution to Si and Fe becomes important towards the center and SN la products have similar 
abundances of Si and Fe at least around M 87, which may reflect dimmer SN la observed in old stellar systems. 
The S abundance is similar to the Si abundance at the center, but has a steeper gradient. This result suggests 
that the S/Si ratio of SN II products is much smaller than the solar ratio. 

Key words. X-rays:ICM — galaxiesdSM — ICM:abundance — individual:M 87 



K>" 1. Introduction 



The intracluster medium (ICM) contains a large amount 
of metals, which are mainly synthesized in early- type 
galaxies (e.g. Arnaud et al. 1992; Renzini et al. 1993). 
Thus, abundances of the metals are tracers of chemical 
evolution of galaxies and clusters of galaxies. 

Based on the Si/Fe ratio observed with ASCA, a dis- 
cussion on contributions from SN la and SN II to the met- 
als has commenced. In a previous nucleosynthesis model 
of SN la, the Fe abundance is much larger than the Si 
abundance in the ejecta of SN la (Nomoto et al. 1984). 
Observations of metal poor Galactic stars indicate that 
average products of SN II have a factor of 2-3 larger abun- 
dance of a-elements than Fe (e.g. Edvardsson et al. 1993; 
Nissen et al. 1994; Thielemann et al. 1996), although this 
ratio may depend on the initial mass function (IMF) of 
stars. From elemental abundance ratios in the ICM of four 
clusters of galaxies observed with ASCA, Mushotzky et al. 
(1996) suggested that the metals in the ICM were mainly 
produced by SN II. Since stars in early-type galaxies in 
galaxy clusters are very old (e.g. Stanford et al. 1995; 
Kodama et al. 1998), this means that most of the met- 
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als in the ICM are produced through star formations at 
high z. Fukazawa et al. (1998) systematically studied 40 
nearby clusters and found that the Si/Fe ratio is lower 
among the low-temperature clusters, which indicates that 
SNe la products are also important among these clusters. 
From the observed radial dependence of the abundances, 
Finoguenov et al. (2000) found that the SNe II ejecta have 
been widely distributed in the ICM. 

In the ICM, around a cD galaxy, the contribution of 
metals from the galaxy becomes important. Fukazawa et 
al. (2000) found a central increment of Fe and Si abun- 
dances around cD galaxies, which is due to SNe la from 
the cD galaxies. The supply of metals by stellar mass loss 
must be also considered (e.g. Matsushita et al. 2000) 

In addition to the Si and Fe abundances, the XMM- 
Newton observatory enables us to obtain a-element abun- 
dances such as for O and Mg, which are not synthe- 
sized by SN la. Bohringer et al. (2001) and Finoguenov 
et al. (2002) analyzed annular spectra of M 87 observed 
by XMM-Newton, and found a flatter abundance gradi- 
ent of O compared to steep gradients of Si, S, Ar, Ca, 
and Fe. The stronger abundance increase of Fe compared 
to that of O indicates an enhanced SN la contribution 
in the central regions (in disagreement with our findings, 
Gastaldello and Molendi (2002) claim a similar gradient 
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for O and Fe, which would not support this conclusion). 
The observed similar abundance gradients of Fe and Si 
and the large Si/0 ratio at the center imply a significant 
contribution of Si by SN la that is a larger Si/Fe abun- 
dance ratio by SN la than obtained by the classical model 
of Nomoto et al. (1984). The larger Si/Fe ratio and the 
implication that the Si/Fe ratio supplied by SN la may 
change with radius in the M 87 halo may indicate a di- 
versity of SN la explosions also reflected in a diversity of 
the light curves as discussed by Finoguenov et al. (2002). 
A similar abundance pattern of O, Si and Fe is observed 
around center of A 496 (Tamura et al. 2001). 

A prerequisite of the abundance determination is a 
precise knowledge of the temperature structure of the 
ICM. Based on the XMM-Newton observation of M 87, 
Matsushita et al. (2002; hereafter Paper I) found that the 
intracluster medium has a single phase structure locally, 
except for the regions associated with radio jets and lobes 
(Bohringer et al. 1995; Belsole et al. 2001), where there 
is an additional ^ 1 keV temperature component. The 
signature of gas cooling below 0.8 keV to zero tempera- 
ture is not observed as expected for a cooling flow (e.g 
Fabian et al. 1984). The fact that the thermal structure 
of the intracluster medium is fairly simple and the plasma 
is almost locally isothermal facilitates the abundance de- 
termination enormously and helps to make the spectral 
modeling almost unique. 

In this paper, based on the detailed study of the tem- 
perature structure, abundances of various elements are 
discussed. In section 2, we summarize the observation and 
data preperation. Sections 3 and 4 deal with the MEKAL 
(Mewe et al. 1995, 1996; Kaastra 1992; Liedahl et al. 1994) 
and the APEC model (Smith et al. 2001) fits to the de- 
projected spectra, and in section 5, abundance ratios are 
determined directly from the line ratios, considering their 
temperature dependences. In section 6, we evaluate an ef- 
fect of resonant line scattering. Section 7 gives a discussion 
of the obtained results, and Section 8 summarizes the pa- 
per. We adopt for the solar abundances the values given 
by Feldman (1992), where the solar Ar and Fe abundances 
relative to H are 4.47x10"^ and 3.24x10"^ in number, 
respectively. These values are different from the "photo- 
spheric" values of Anders & Grevesse (1989), where the so- 
lar Ar and Fe abundances are 3.63x10"^ and 4.68 xlO^'^, 
respectively. The solar abundances of the other elements 
are consistent with each other. Unless otherwise specified, 
we use 90% confidence error regions. 

2. Observation 

M 87 was observed with XMM-Newton on June 19th, 
2000. The effective exposures of the EPN and the EMOS 
are 30ks and 40ks, respectively. The details of the anal- 
ysis of background subtraction, vignetting correction and 
deprojection technique are described in Paper I. When 
accumulating spectra, we used a spatial filter, excluding 
those regions where the brightness is larger by 15 % than 
the azimuthally averaged value in order to excise the soft 



emission around the radio structures (Figure 1 and 6 in 
Paper I), although it is not fully excluded within 2', due 
to its complicated structure and the limited spatial reso- 
lution of the XMM telescope. 

For the EPN spectra, we employ the response ma- 
trix corresponding to the average distance from the read- 
out for each accumulated region, epn_fs20_sYi_thin.rmf 
from March 2001. Here, i reflects the distance from the 
readout-node, which affects the energy resolution. For the 
EMOS data, we used ml_thinlv9ql9t5r5_all_15.rsp from 
June 2001. The spectral analysis uses the XSPEC.vll.l 
package. We fitted the EPN and EMOS data in the spec- 
tral range of 0.5 to 10 keV. For outside 8', the energy range 
between 7.5 and 8.5 keV of the EPN spectra is ignored, 
because of strong emission lines induced by particle events 
(Freyberg et al. 2002). 

3. Spectral fit with the MEKAL model 

We have fitted the deprojected spectra with a single tem- 
perature MEKAL model with photoelectric absorption, in 
the same way as in Paper I. For the outermost region, we 
fitted the projected (annular) spectrum. For the spectra 
within 1', a power-law component is added for the central 
active galactic nucleus. We fixed index of the power-law 
component to the best-fit value obtained from the spec- 
trum within 0.12' (Paper I) and normalized it using the 
point spread function of the XMM telescope. Abundances 
of C and N are fixed to be 1 solar and those of other ele- 
ments are determined separately. Within 2', we also fitted 
the spectra with a two temperature MEKAL model, where 
the abundances of each element of the two components are 
assumed to have the same values, since within this radius, 
there remains a small amount of the cooler component 
with a nearly constant temperature of 1 keV associated 
with the radio lobes (Belsole et al. 2001; Paper I). As 
discussed in Paper I, the two temperature model should 
refiect the actual temperature structure of the ICM, al- 
though within 0.5', the XMM spatial resolution is not 
enough to resolve temperature components of the com- 
plicated structure of the ICM. 

The results are summarized in Table 1 and Figures 1- 

4. Representative spectra fitted with the MEKAL model 
are shown in Figure 2 and 7 of Paper I and Figure 0. There 
remain small discrepancies between the data and the best 
fit model. A residual structure exists at 0.8-1 keV of the 
deprojected spectra of the EPN of R >2', when fitted with 
the single temperature MEKAL model, while the EMOS 
data have no such residual (Figure 2 of Paper I). Here, R is 
the 3-dimensional radius. This problem will be discussed 
in the next paragraph. The Fe-L/Mg-K structure around 
1.3 keV is not fitted well (Figure |l]). Therefore, we do not 
show Mg abundances here and will discuss this in more 
detail below (Sec. 4 and 5.5-5.6). We also do not show 
Ne and Ni abundances, because their line strengths are 
much smaller than those of the Fe-L features at the same 
wavelength and a small uncertainty in the Fe-L complex 
gives large uncertainties in abundances of these elements 
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Table 1. Results of spectrum fitting of the deprojected spectra with the single temperature and the two temperature 
MEKAL model. For the outermost region, projected spectra are used. The regions with the enhanced X-ray emission 
associated with the radio structures are excluded throughout this paper. 

R kTl kT2 iVn O Si S Ar Ca Fe xVm" 

(arcmin) (keV) (keV) (lo^^cm^^) (solar) (solar) (solar) (solar) (solar) (solar) 

The EMOS spectra fitted with a IT MEKAL model 

0.00- 0.35 ETS O 002 037 037 Tol L25 O20 435/177 
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channel energy (keV) 

Fig. 1. The data-to-model ratios of the deprojected spec- 
tra of the EMOS within R=2-S' fitted with the single tem- 
perature MEKAL model (black crosses) and version 1.1 of 
the APEC model (gray diamonds). 



(Masai 1987; Matsushita et al. 1997; 2000). There are also 
small discrepancies in the continuum level around Ka lines 
of Si, S and Ar, which will be discussed in Sec. 5. 

The derived abundances are mostly consistent between 
the EPN and the EMOS, although the EPN abundances of 
O, Si, S and Fe are systematically smaller than the EMOS 



results by 10-30%, when fitted with a single temperature 
MEKAL model. As will be shown in Sec. 5, the strengths 
of Ka fines of H-fike O, Si, S and He-like (-1- Li fike) Fe 
observed by the EPN and EMOS agree within several per- 
cent when considering the difference in the normalization 
between the two detectors. Therefore, these abundance 
discrepancies should be caused by the discrepancy at 0.8- 
1 keV, which may reflect uncertainties in an instrumental 
low energy tail of the strong Fe-L fines observed in EPN 
spectra, as discussed in Paper I. In addition, the EPN is 
characterized by a dependence of the energy resolution 
on the position on the detector. As a result, the observed 
Fe abundance of the EPN may be artificially decreased 
by these problems and the change of the continuum level 
also affects the O abundance. Those of Si and S also cou- 
ple with abundances of other elements such as O and Fe 
due to the effect of free-bound emission and the increment 
of the strength of the Fe-L emission around the Ka line of 
He-like Si. When we fix the O, Ne, Mg, Fe and Ni abun- 
dances to the best fit values derived from the EMOS spec- 
tra of the whole energy band and fitted the EPN spectra 
above 1.7 keV, we recover consistent Si and S abundances 
between the two detectors (Figure ||) . 

The derived abundances of Si, S, Ar, Ca, and Fe from 
deprojected spectra using the single temperature MEKAL 
model fit show strong negative gradients at i? > 1', and 
drop sharply within this radius (Table 1, Figure ||) as al- 
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Fig. 2. Abundance profiles of O, Si, S, Ar, Ca and Fe using the single temperature (black) and the two temperature 
(gray) MEKAL model from the deprojected data of the EMOS (diamonds) and the EPN (closed circles). The open 
circles are the Si and S abundances derived from the EPN spectra above 1.7 keV, when O, Ne, Mg, Fe and Ni 
abundances are fixed to be the best fit values derived from the EMOS data of the whole energy band. The dashed 
lines correspond to the best fit relation of the deprojected Si profile of the EMOS. 



ready seen in the projected spectra (Bohringer et al. 2001; 
Finoguenov et al. 2001; Gastaldello & Molendi 2002). But 
now the absolute values are slightly larger than those from 
projected spectra. In contrast, the O abundance is nearly 
constant at i? > 1'. 



The two temperature MEKAL fit on the deprojected 
spectra gives significantly larger abundances than the sin- 
gle temperature MEKAL fit (Figure |[]|, Table 1) and the 
central abundance drops seen in the single MEKAL fit 
become very week in the results of the two temperature 
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MEKAL fit, although contributions of the cooler compo- 
nent are small as shown in Paper I. For example, at i?=0.5- 
1', 10 percent of the cooler component in units of emission 
measure changes the Fe abundances by a factor of 2. The 
reason for this is the change of the temperature of the hot 
component, when the cold component is added. For a fixed 
Fe line feature this leads to a higher abundance required 
in the fit. In addition, abundances of O, Si and S also in- 
crease by a factor of 2. The increment of the O abundance 
is due to the change of the temperature structure. The 
Si and S abundances increase due to an increment of the 
free-bound emission and the strength of Fe-L lines around 
the He-like Si line. Therefore, considering the remaining 
1 keV temperature component associated with the radio 
structures within 2', abundances become nearly constant 
within 2'. 

We have also tried a three temperature MEKAL model 
for spectra within 2'. Although the derived values have 
slightly improved, the derived abundances and their er- 
rors have nearly the same values as those from the two 
temperature MEKAL model fit. 

Figure || shows the confidence contours for the Si, S 
and Fe abundances, derived from the deprojected spec- 
tra. For the outermost region, the results from the pro- 
jected spectra are plotted. The elliptical shape of the con- 
fidence contours indicates that abundance ratios are better 
determined than the abundances themselves. The differ- 
ences in the Si/Fe and S/Fe ratios between the EMOS and 
EPN are considerably smaller than the abundance differ- 
ences. Furthermore, the two temperature MEKAL model 
fit gives almost the same values of the Si/Fe and S/Fe ratio 
with the single temperature MEKAL model fit, although 
absolute values of Si, S, and Fe are largely changed. This 
means that the abundance ratios do not strongly depend 
on the temperature structure. As a result, the Si/Fe ra- 
tio is determined to be nearly constant at 1.1± 0.1 solar, 
with no radial gradient. In contrast, we find that the S 
abundance has a steeper gradient than Fe abundance. 

The abundances of Si and Fe are nearly independent 
to the O abundance (Figure However, at the center, 
the two temperature MEKAL model gives the similar val- 
ues of the 0/Si and 0/Fe ratio as the single temperature 
MEKAL model. The 0/Si and 0/Fe ratios are consistent 
between the two detectors, although the abundances of O, 
Si, and Fe between the EMOS and EPN differ by 10-30%. 
Both the Si/0 and Fe/0 ratios increase by a factor of 1.5 
with radius. This result is in disagreement with the claim 
by Gastaldello & Molendi (2002), where they derived the 
constant Fe/0 ratio. The background subtraction should 
not be a problem when determining O, Fe, and Si abun- 
dances, since even at the radius of 10', the background 
is only a few percent between 0.5 to 2 keV. In addition, 
the background was subtracted according to the detailed 
study by Katayama et al. (2002). A more severe problem 
should be the uncertainties in the response matrix, since 
the O abundance was inconsistent between the EPN and 
the EMOS by a factor of 2 using an old response matrix 
(Bohringer et al. 2001), although the discrepancy becomes 
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Fig. 3. Confidence contours (90%) for the Si, S vs. Fe 
abundance of the deprojected data of the EMOS (thick 
lines) and the EPN (thin lines) fitted with the single tem- 
perature MEKAL model (solid lines) and the two temper- 
ature MEKAL model (dotted lines). For the outermost 
region, the projected spectra are used. The solid lines in- 
dicate the condition of the solar-abundance ratio. The 
dashed lines represent the Si/Fe ratio to be 0.8 and 1.2 
solar and the dotted line represents the S/Fe ratio to be 
0.5 solar. 



smaller now as shown in this section. Using an old version 
of the response matrix and fitting the projected spectra in 
the same way in Gastaldello & Molendi (2002), we achieve 
the same results. 

In summary, from the deprojected spectra using the 
MEKAL model, the abundances of Si, S, Ar, Ca, and Fe 
have steep gradients at i? > 2' and become nearly con- 
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Fig. 4. 90 % confidence contours for the Si, Fe vs. O. The 
meanings of the symbols are the same as in Fig. 3. The 
sohd lines and dot-dashed lines represent the Si or Fe to 
O ratio to be 1 and 2 solar, respectively. 

stant within the radius. The Si/Fe ratio is almost con- 
stant within the field of view, while the S/Fe ratio has a 
negative, and the 0/Fe ratio has a positive gradient. 

4. Spectral fits with the APEC model 

A new plasma code, APEC (Smith et al. 2001), is now 
available in XSPEC_vll. Within the energy resolution of 
the CCD detectors, the dominant change from version 1.0 
to the APEC code in the XSPEC version 11.0.1 is not 
in the Fe-L lines, but in the Ka lines. Using version 1.0 
of the APEC code, strengths of the Ka lines of H-like 
ions, which are the most fundamental lines, decrease by 
~ 50 % at 2 keV, and their temperature dependences also 



change. The strength of the 6.7 keV Fe-K line decreases, 
by a factor of 2 at 1 keV and 1.3 at 4 keV. The changes 
of the Ka lines of He-like ions of Si, S, Ar and Ca are 
smaller. As a result, any single temperature APEC model 
cannot fit the deprojected spectra of M 87, especially for 
the ratio of the Fe-L and the Fe-K as shown in Paper I, 
although a single temperature MEKAL model can fit each 
deprojected spectrum at i? > 2'. 

However, the strengths of Ka lines of version 1.1 of the 
APEC code in version 11.1 of the XSPEC became more 
consistent with the MEKAL ones. Within the temperature 
range of the ICM of M 87, the largest differences except 
for the Fe-L structure are 20% changes in the strengths of 
the Ka line of H-like O and the Fe line at 6.7 keV. 

Using version 1.1 of the APEC code, we fitted the de- 
projected spectra of the EMOS. The results are summa- 
rized in Table 2 and Figure ^ In contrast to the worse 

values obtained from version 1.0 of the APEC code, 
version 1.1 gives similar values compared to those ob- 
tained by the MEKAL model. For the Fe-L/Mg-K struc- 
ture around ~ 1.4 keV, the APEC (vl.l) model gives a 
better fit than the MEKAL model (Figure |l]). The derived 
temperatures and hydrogen column densities are consis- 
tent with those from the MEKAL model. The contribu- 
tion of the lower temperature component at i?=l-2' has 
sHghtly decreased from 4% by the MEKAL code to 2% 
by the APEC code, although the contribution within 1' 
is consistent with each other. The few percent of the soft 
component using the MEKAL model at i?=l-2' may be 
an artifact due to the uncertainties in the Fe-L complex. 

Figure |^ compares the abundances derived from the 
APEC and the MEKAL model. The large abundance 
changes seen in APEC version 1.0 (Paper I) are not de- 
rived using the new APEC code, version 1.1. For Si, S, and 
Fe, we find consistent results within several percent. The 
O abundances from the APEC code are systematically 
larger than those from the MEKAL code by 20%, which 
is close to the difference of the line strength of the Ka line 
of H-hke O between the MEKAL and APEC model. 

In contrast, the derived Mg abundances differ by a 
factor of 2 between the two codes. Since the APEC model 
can fit the Mg-K/Fe-L structure at ^^1.4 keV quite well, 
the APEC code may be better to derive a Mg abundance. 
There is another problem due to the strong instrumental 
line of the EMOS, which is located at a similar energy to 
the Mg line. At r=8-13.5', the strength of the instrumen- 
tal line is a factor of 10 larger than the Mg line of the 
ICM, and the uncertainty affects the spectrum at i?=4-8' 
through the deprojection. Here, r is the projected radius 
from the center. Therefore, the Mg abundances are shown 
only within 4'. The derived Mg abundance is 1 solar 
at i? < 2' and shows a radial gradient outside this radius 
(Table 2). 

The observed values of the Ne abundance derived from 
the APEC model are also a factor of 1.3 larger than those 
derived from the MEKAL model. For the Ni abundance, 
the APEC code gives 60% of the values from the MEKAL 
code at i? > 2', but within R < 2' , both codes give similar 
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Table 2. The results of the spectral fits of the deprojected spectra of the EMOS using the APEC (vl.l) code 



R 


kTl 


kT2 


Ah 





Mg 


Si 


s 


Fe 


2 / a 

X 


(arcmin) 


(keV) 


(keV) 


(lO^^cm^^) 


(solar) 


(solar) 


(solar) 


(solar) 


(solar) 




a IT APEC model 


0.00- 0.35 


1.15 




2.5 


0.06 


0.0 


0.29 


0.41 


0.15 


253/175 


0.35- 1.00 


1.48 




2.4 


0.97 


1.17 


1.80 


1.79 


1.37 


583/175 


1.00- 2.00 
2.00- 4.00 
4.00- 8.00 
8.00-11.30 
11.30-13.50 


1 50+0-01 
^•'-'J^-o.oi 

9 90+0.02 
^•^°-0.03 
9 ^0+0.03 

9 7/1+0.05 
^- '^-0.07 




2.811:? 

9 1+0.3 
^^■^-0.3 

2 0+°-^ 

^•'J-0.4 
9 9+0.7 
^•^-0.7 

0.8«-j 


56+° " 

U.OD_o og 
^■3*-0.09 

0.52t°:°^ 
0.46t°:- 

0.38t°:J^ 


0.87t°;i^ 
0.67t°:i? 


1 47+0-10 
1 iq+0-07 

^•^''-0.06 

0.871°°^ 
0.65«:;« 
0.62t°°^ 


i.28t°:- 

0.68«:°^ 
iq+0-09 

"•^^-0.08 


1 -r p+0.05 

0.78lH^2 

64+° °*' 

'-'•"^-0.05 

5S+° °=^ 


126/175 
181/175 
210/175 
248/175 
305/175 


a 2T APEC model 


0.00- 0.35 
0.35- 1.00 
1.00- 2.00 


1 7S+0.06 

1 fiS+°°l 
i-DO_0.oi 


'J-'°-0.02 
J^-'JJ_0.02 

1.00 (fix) 


2 4+^-^ 

2.o«:« 

2.5«-i 


0.26t°:f6 

o.8ot°:« 

0.64t°:^g 


0.75tr3( 

i-ioi°:i^ 
i.ooi°:f. 


1 14+0-9!) 

1 70+O-O9 
^•'^-0.16 

J^-DJ^-0.23 


1 -10+0.03 
i-i»_0.42 

1 61+°" 

1 00+0. 26 
i.OO_Q 22 


96+°'^° 

1 58+0-06 

1 26+°" 


138/174 
190/174 
125/174 



° Degrees of freedom 
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Fig. 5. The ratios of the best fit abundances derived from 
the APEC model and MEKAL model. The two tempera- 
ture model is used at i? < 2' and the single temperature 
model is used at i? > 2'. 

values. In contrast to the Mg lines, the Ni-L lines and the 
Ne-K lines cannot be distinguished in the spectra of M 87. 
Therefore, we do not present Ni and Ne abundances in this 
paper. 

In summary, the new APEC code provides a better fit 
to the Fe-L structure, and the derived abundances of O, 
Si, and Fe are mostly consistent with those obtained from 
the MEKAL code. The Mg abundances are derived to be 
~ 1 solar at i? < 2', since the Fe-L/Mg-K structure is well 
fitted with the APEC code. 

5. Abundance determination using line ratios 

In this section, we try to obtain abundance ratios of el- 
ements directly from ratios of emission lines. When cal- 
culating the deprojected spectra, we make an assumption 
that spectra beyond the field of view are the same as the 
outermost spectrum, which is not valid since the temper- 
ature starts to decrease beyond the field of view of the 



detector (Shibata et al. 2001). Therefore, especially in the 
outer regions, there may be a bias due to this assump- 
tion, since emission lines have a stronger dependence on 
temperature than the continuum. In addition, abundances 
from deprojected spectra have a larger uncertainty due to 
poorer statistics than those derived from projected data. 
In order to derive abundances more accurately, we need 
projected spectra. 

However, the projected spectra consist of multi- 
temperature components. As a result, some line ratios, 
such as the ratio of Ka lines of He- and H-like S, of the 
projected spectra cannot be fitted with a single temper- 
ature MEKAL model (Figure ||). We must be careful in 
using a multi-temperature model, e.g. a two temperature 
model, since emission lines and a continuum spectrum are 
different functions of temperature. In addition, the contin- 
uum spectra between S, Ar and Ca lines show small dis- 
crepancies between the data and the single temperature 
model. These small discrepancies sometimes give larger 
uncertainties in the strengths of faint lines such as Ka 
lines of He-like Ar and Ca. For example, when fitting with 
a single temperature MEKAL model with the projected 
spectrum at 2-4', temperature and Ar abundance are de- 
termined as 2.2 keV and 0.7 solar, respectively. However, 
when we restrict the energy range within 2.8-5.0 keV and 
fix the abundances of O, Ne, Mg, Si, S and Fe to the best 
fit values obtained from the whole energy band fit, the 
temperature becomes 2.4 keV and Ar abundance is de- 
termined as 1.1 solar. Here, the Ar abundance is changed 
by 60%, although the best fit temperatures differ by only 
10%. This small discrepancy between the S and Ar lines is 
also seen in the deprojected spectra, although their errors 
are larger. 

Therefore, in this section, we determine strengths of 
emission lines, and then determine abundance ratios, 
considering the temperature dependence of line ratios. 
We have fitted the projected and deprojected spectra 
within an energy band around a given line with thermal 
bremsstrahlung and gaussians. Obtaining strengths of Ka 
lines of He-like S and Ar, contributions of K/3 lines of H- 



8 Kyoko Matsushita et al.: XMM observation of M 87 



Table 3. Ratios of Ka lines and abundance ratios (errors correspond to 68%). Within 2', the effect of the lower 
temperature component is considered. 



r 




H-O/H-Si 


0/Si (solar) 


(H-S±He-S)/H-Si 


S/Si (solar) 


He-Ar/H-Si 


Ar/Si (solar) 


He-Ca/H-S 


Ca/S (solar) 


Annular spectra of the EMOS 


0.0- 


0.5 


2.26 ±0.29 


0.43 ±0.06 


1.08 ±0.06 


1.14 ±0.06 


0.16 ±0.03 




0.20 ±0.04 




0.5- 


1.0 


2.16 ±0.15 


0.48 ± 0.03 


1.02 ± 0.04 


1.05 ±0.05 


0.14 ±0.02 


0.79 ±0.11 


0.23 ± 0.03 


1.29 ±0.16 


1.0- 


2.0 


1.90 ±0.11 


0.45 ± 0.03 


0.94 ± 0.03 


0.96 ± 0.03 


0.18 ±0.03 


0.87 ±0.12 


0.28 ± 0.03 


1.53 ±0.15 


2.0- 


4.0 


1.84 ±0.13 


0.46 ± 0.03 


0.93 ±0.03 


0.97 ±0.03 


0.17 ±0.01 


0.78 ± 0.06 


0.23 ± 0.02 


1.16 ±0.10 


4.0- 


8.0 


2.38 ±0.10 


0.59 ±0.02 


0.82 ±0.03 


0.87 ±0.03 


0.14 ±0.02 


0.67 ±0.09 


0.24 ± 0.02 


1.18 ±0.11 


8.0-13.5 


2.45 ±0.20 


0.61 ±0.05 


0.62 ±0.05 


0.67 ±0.05 


0.14 ±0.03 


0.71 ±0.16 


0.32 ± 0.06 


1.56 ±0.29 



Deprojected spectra of the EMOS 



0.0- 


0.5 


2.08 ±0.42 


0.35 ±0.07 


1.09 ±0.17 


1.16 ±0.18 


0.18 ±0.06 




0.18 ±0.08 




0.5- 


2.0 


2.06 ±0.23 


0.45 ±0.05 


0.96 ± 0.04 


1.00 ±0.04 


0.16 ±0.02 


0.93 ±0.11 


0.29 ± 0.04 


1.60 ±0.21 


2.0- 


5.6 


1.77 ±0.12 


0.44 ± 0.03 


0.90 ± 0.04 


0.93 ± 0.04 


0.18 ±0.02 


0.85 ±0.10 


0.24 ± 0.03 


1.24 ±0.16 


5.6-11.3 


2.74 ± 0.31 


0.68 ± 0.08 


0.75 ± 0.07 


0.80 ± 0.08 


0.14 ±0.03 


0.67 ±0.15 


0.25 ± 0.05 


1.26 ±0.26 


Annular spectra of the EPN 


0.0- 


0.5 


2.00 ±0.33 


0.38 ±0.06 


1.11 ±0.11 


1.17 ±0.11 


0.30 ± 0.05 




0.24 ±0.09 




0.5- 


1.0 


1.87 ±0.20 


0.41 ±0.04 


1.04 ± 0.06 


1.07 ±0.06 


0.15 ±0.03 


0.81 ±0.13 


0.24 ± 0.06 


1.36 ±0.31 


1.0- 


2.0 


1.90 ±0.18 


0.44 ± 0.07 


0.92 ±0.05 


0.94 ±0.05 


0.18 ±0.02 


0.88 ±0.10 


0.28 ± 0.04 


1.50 ±0.22 


2.0- 


4.0 


1.87 ±0.23 


0.47 ±0.06 


0.93 ±0.05 


0.97 ± 0.05 


0.17 ±0.02 


0.81 ± 0.11 


0.35 ± 0.05 


1.78 ±0.25 


4.0- 


8.0 


2.09 ±0.26 


0.52 ±0.06 


0.76 ± 0.05 


0.82 ± 0.05 


0.08 ± 0.02 


0.41 ±0.09 


0.34 ± 0.09 


1.69 ±0.44 


8.0-13.5 


2.70 ± 0.43 


0.67 ±0.11 


0.67 ±0.08 


0.73 ± 0.09 


0.20 ±0.05 


1.02 ±0.23 


0.44 ± 0.08 


2.14 ±0.39 



m 



a 




2.5 3 



channel energy (keV) 

Fig. 6. The projected EMOS spectrum of r=8.0-13.5' 
(crosses) fitted with a single MEKAL model (gray) and 
bremstrallung and gaussians (black). 



like Si and S were subtracted using the strengths of Ka 
lines. Here, we used the K/3 to Ka ratio of the MEKAL 
code. Strengths of H-like Ka lines of O were obtained from 
spectra within the energy range of 0.55-1 keV fitted with 
a gaussian and two APEC models with zero O abundance. 
Those of Si were also obtained from fitting with gaussians 
and the two temperature APEC model with zero Si abun- 
dance, using the 1.1-2.2 keV energy band. We selected 
line ratios whose temperature dependence is small, for ac- 
curate measurements of abundance ratios. The results are 
summarized in Table 3. 



5.1. vs. Si 

The abundance ratio of 0/Si is better obtained from the 
ratio of the line strengths of Ka lines of H-like ions. 
Although the O and Si abundances obtained from the 
EMOS and the EPN using the same single temperature 
model have small discrepancies of 20~ 30%, the line ra- 
tios of the two detectors agree well with each other (Figure 
I). 

The observed radial profile of the line ratio has a mini- 
mum at ~2' (Figure]^). Figure ^also shows the line ratios 
plotted against the best fit MEKAL temperatures from 
the whole energy band. When the abundance ratio of 0/Si 
is constant, the line ratio is almost constant above 1.7 keV 
and it starts to increase sharply below the lowest temper- 
ature. Comparing the observed minimum value of ~ 2 
with the minimum value of ^ 4 corresponding to the solar 
abundance ratio, no temperature distribution can repro- 
duce an 0/Si ratio larger than 0.5 solar at the radius of 
~2'. 

Considering the temperature structure derived in 
Paper I, we can better constrain the 0/Si ratio. The in- 
dependence of the line ratio from the temperature above 
1.7 keV means that the abundance ratio can be directly 
obtained dividing the line ratio by that of the solar abun- 
dance ratio when there is no temperature component be- 
low 1.7 keV. Thus, beyond 2', the gradient of the line ratio 
reflects a change of the 0/Si ratio, since outside 2', there 
is no temperature component below 1.7 keV (Paper I). 
From the line ratio at r= 2-4', the 0/Si is derived to be 
0.46 solar, and beyond r > 8', it increases to 0.6-0.7 solar 
(Figure 0). 

Within 2', the presence of the additional ~ 1 keV 
temperature component must be considered. The derived 
0/Si ratio assuming that the temperature is larger than 
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1.7 keV should be the maximum value of the 0/Si ratio 
which is ^ 0.5 solar within 2' (triangles in the bottom 
panel of Figure]^). We also derive the 0/Si ratio from the 
line ratio considering the fraction of the low temperature 
component, which is derived from the spectral fit using the 
best fit two temperature APEC model. Then, the 0/Si ra- 
tio is derived to be 0.45 solar at 1-2' and 0.4 solar within 
0.5'. Since the line ratio is a steep function when the tem- 
perature is around 1 keV, we have also checked a case 
when the lower temperatures is 0.6 keV, which is 0.2 keV 
smaller than the central temperature. For the innermost 
region, the derived 0/Si ratio with the temperature of 0.6 
keV is consistent within 10 percent, since the temperature 
of the softer component becomes lower, the fraction of the 
component also becomes lower which is derived from the 
spectral fitting from Fe-L. 

For comparison, we also plotted the 0/Si ratio derived 
from the deprojected spectra through the spectral fits and 
the line ratio. Although the 0/Si ratios derived from the 
line ratios are systematically larger by ~ 20% than those 
derived from spectral fit with the MEKAL model, the gra- 
dient of the 0/Si ratio is consistent between the two, and 
the 0/Si ratio changes by a factor of 1.4 from the center 
to 10'. 



5.2. S vs. Si 

Although the ratios of Ka lines of He-like or H-like S to 
that of H-like Si are both steep functions of temperature, 
the ratio of the sum of the two Ka lines of S to the H-like 
Si line is nearly constant within 10% between 0.9 to 5 keV 
when the S/Si ratio is constant (Figure ||). The important 
point is that we can derive the abundance ratio almost in- 
dependently from the temperature structure of the ICM. 
Considering that there is no temperature component be- 
low 1 keV, as in the case of the 0/Si ratio, even for the 
projected data, the S/Si ratio can be directly calculated 
from the line ratio. As in the previous subsection, the ef- 
fect of the lower temperature component is estimated to 
be less than a few percent. 

The observed profile of the line ratio shows a negative 
gradient and at r=10', it is about half of the central value. 
This gradient reflects the change of the S/Si ratio, because 
of the small dependence on temperature. Converted to the 
abundance ratio, the S/Si ratio is 1.0 solar within 2', and 
drops to 0.67 solar at ^ 10' (Figure ||). These results are 
systematically larger than those derived from the spectral 
fits. Especially, at r > 8', the spectral fits on the projected 
spectra give the S/Si ratio to be 0.3-0.5 solar (Table 1, 
Figure 2), while the line ratio gives the value of 0.6-0.7 
solar. For the projected data, the S abundance from the 
spectral fit should not be correct since the single temper- 
ature model cannot fit the two Ka lines simultaneously as 
shown in Figure 6. 
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Fig. 7. The line ratios of Ka lines of H-like O and Si of the 
EMOS (black) and the EPN (gray) , plotted against the ra- 
dius (upper panel) and the best fit temperature fitted with 
the single MEKAL model (middle panel). Solid and dot- 
ted data correspond to projected and deprojected data, 
respectively. Dashed lines in the middle panel correspond 
to constant abundance ratios in unit of the solar ratio. 
The derived abundance ratios are plotted in the bottom 
panel. Within 2 arcmin, contributions of the lower tem- 
perature component are considered. Triangles are the best 
fit values of the abundance ratio, assuming the tempera- 
ture is larger than 1.7 keV. The diamonds correspond to 
the abundance ratios derived from the deprojected EMOS 
snectra usinff the MEKAL model fit. Errors renresent the 
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Fig. 8. The line ratios of Ka lines of the sum of Hc-likc 
and H-hke S to H-like Si of the EMOS (black) and the EPN 
(gray), plotted against radius (top panel) and temperature 
(middle panel), and the radial profile of the S/Si ratio 
(bottom panel). The meaning of the symbols are the same 
as in Figure 7. Errors represent the 68% confidence level. 
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Fig. 9. The line ratios of Ka lines of He-like Ar to H-like 
Si of the EMOS (black) and the EPN (gray) (top panel) 
and the Ar/Si ratio (bottom panel). The meaning of the 
symbols are the same as in Figure 7. Errors represent the 
68% confidence level. 

Ar/Si, the change of the line ratio is less than 10 % (Figure 
|). In Paper I, we found that at i? > 0.5', the ICM tem- 
perature at a given radius is dominated by a single tem- 
perature component whose temperature is larger than 1.7 
keV. From Figure 7 in Paper I, the contributions from the 
1 keV component to Ar and H-like Si lines are negligible. 
Therefore, as in case of the S to Si abundance ratio, we can 
obtain the abundance ratio of Ar/Si from the line ratio. 
The results indicate that the Ar/Si is 0.7-0.8 solar. 

The derived Ar/Si ratios are systematically larger 
than the results of the spectral fit on the projected data 
(Finoguenov et al. 2002) by 30% within 3' and by a factor 
of 2 at r > 8'. This large discrepancy is caused by the 
failure to fit the continuum between the S and Ar lines 
illustrated in Figure 6. 



5.3. Ar vs. Si 

Although the line ratio of the Ka lines of He-like Ar to 
H-like Si is a steeper function of temperature than S/Si 
or 0/Si ratio, within 1.7 to 2.8 keV, at a given ratio of 



5.4. Ca vs. S 

The ratio of the Ka line of the He- like Ca to that of the H- 
like S also show a small temperature dependency (Figure 
|l0| ) . Because these lines are dominated by a hotter temper- 
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Fig. 10. The line ratios of Ka lines of He-like Ca to H-like 
S of the EMOS (black) and the EPN (gray). The mean- 
ing of the symbols are the same as in Figure 7. Errors 
represent the 68% confidence level. 



ature component even around the center, we can calculate 
the Ca/S ratio in the same way. The Ca/S ratio is ^ 1.5 
solar within the whole region. 

5.5. Mg vs. O 

The energy of the Ka line of H-like Mg is slightly shifted 
from the peak of the Fe-L emission at ^ 1.49 keV (Figure 
|ll|).Thus, the Ka line can be distinguished within the en- 
ergy resolution of the CCDs. The problem with the spec- 
tral fitting is the modeling of the Fe-L structure around 
the Mg line. Due to the problem of the strong instrumen- 
tal line of the EMOS, we do not present the results on the 
Mg-K line for r > 8'. 

The strengths of H-like and He-like Ka lines of Mg 
were obtained from spectra within 1.1-1.6 keV, fitted with 
gaussians and a two temperature MEKAL or APEC model 
with zero Mg abundance. The results are summarized in 
Table 4. 

Figure |l^ shows the ratios of the Ka line of the H-like 
Mg to that of H-like O. The strengths of the Mg line ob- 
tained using the APEC Fe-L model are ~ 30% larger than 
those derived using the MEKAL model. Since the APEC 
fit to the Fe-L/Mg-K feature of 1.4 keV is better than the 
MEKAL model fit, the APEC code may be better to de- 
scribe the Fe-L feature, although the discrepancy between 
the data and model of the MEKAL plus gaussians is only 
several percent. 

The projected EMOS spectra in the energy range of 
1.35 to 1.6 keV are also fitted with a bremsstrahlung model 
plus two gaussians at 1.47 keV and 1.49 keV, since the 
Ka line of H-like Mg at 1.47 keV and a Fe-L peak at 1.49 
can be distinguished within the EMOS energy resolution 
(Figure ^l]). Freeing the strength of 1.49 keV peak of the 
Fe-L enable us to constrain the Mg line strength. The 
derived strengths of the Mg line are consistent with those 




1.3 1.4 1.5 1.6 

channel energy (keV) 

Fig. 11. The projected spectra of r—2-A' at the energy 
band of Mg of the EMOS and the EPN fitted with the two 
temperature APEC model with zero Mg abundance (black 
dotted lines) and gaussians. Those using the MEKAL 
model are shown as gray lines. 

Table 4. Line ratio of Ka of H-like Mg, O, and Si and 
their abundance ratios, when the Fe-L is modeled with the 
APEC (errors correspond to 68%) 



r 




H-Mg/H-O 


Mg/0 (solar) 


H-Mg/H-Si 


Mg/Si (solar) 


Annular spectra of the EMOS 


0.0- 


0.5 


0.23 ±0.04 


1.46 ±0.23 


0.53 ± 0.05 


0.63 ± 0.07 


0.5- 


1.0 


0.22 ±0.02 


1.33 ±0.11 


0.48 ± 0.03 


0.64 ± 0.04 


1.0- 


2.0 


0.21 ±0.02 


1.26 ±0.14 


0.40 ± 0.04 


0.57 ±0.06 


2.0- 


4.0 


0.22 ±0.02 


1.30 ±0.14 


0.40 ± 0.03 


0.59 ± 0.05 


4.0- 


8.0 


0.19 ±0.02 


1.16 ±0.10 


0.46 ± 0.04 


0.69 ± 0.06 






Deprojected spectra 


of the EMOS 




0.0- 


0.5 


0.27 ±0.07 


1.73 ±0.48 


0.56 ±0.13 


0.61 ±0.14 


0.5- 


2.0 


0.20 ±0.03 


1.19 ±0.17 


0.41 ± 0.04 


0.54 ± 0.05 


2.0- 


4.0 


0.20 ±0.05 


1.21 ±0.28 


0.42 ± 0.08 


0.66 ±0.10 


Annular spectra of the EPN 


0.0- 


0.5 


0.31 ±0.06 


1.90 ±0.36 


0.61 ±0.08 


0.72 ±0.10 


0.5- 


1.0 


0.22 ±0.04 


1.31 ±0.22 


0.41 ±0.06 


0.54 ± 0.07 


1.0- 


2.0 


0.19 ±0.03 


1.14 ±0.20 


0.36 ± 0.06 


0.52 ± 0.08 


2.0- 


4.0 


0.23 ±0.05 


1.41 ±0.31 


0.44 ± 0.08 


0.65 ±0.13 


4.0- 


8.0 


0.24 ±0.06 


1.45 ±0.38 


0.50 ±0.12 


0.75 ±0.18 


8.0-13.5 


0.15 ±0.04 


0.94 ±0.26 


0.41 ±0.10 


0.63 ±0.15 



derived using the APEC code for the modeling of the Fe- 
L. This result supports the Mg line strengths from the 
APEC code than those from the MEKAL code. 

As in Sec. 5.1-5.4, we have plotted the line ratio from 
the APEC code against the best fit MEKAL temperature 
(Figure [l2|) . The constant Mg/0 ratio gives a constant 
line ratio within 10% above 1 keV. As in Sec. 5.1, we 
have estimated the effect of the temperature component 
below 1 keV, which is less than 10% for the innermost 
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Fig. 12. (Top panel) The ratio of strength of Ka hnes 
of H-Uke Mg and O of the projected spectra of the 
EMOS (black) and EPN (gray). The Fe-L is modeled with 
the APEC (diamonds) and the MEKAL (crosses) and 
bremsstrahlung and gaussians (bold dashed lines). The 
results of the deprojected spectra of the EMOS, where 
the Fe-L lines are modeled with the APEC are shown as 
dotted diamonds. (Bottom panel) The line ratios when the 
Fe-L lines are modeled with the APEC are plotted against 
the best fit temperature fitted with the single MEKAL 
model. Dashed lines correspond to constant abundance 
ratios. Errors represent the 68% confidence level. 



region. The derived Mg/0 ratios are about 1.2-1.3 solar 
with no radial gradient at r > 0.5', when we adopt the 
Fe-L modeling for the APEC (Figure 13). Within 0.5', a 
higher Mg/0 ratio is allowed, since adding a temperature 
component blow 1 keV, the Mg/0 ratio increases. 

The O and Mg abundances obtained from the RGS 
spectrum through a spectral fit are 0.49±0.04 solar and 
0.9±0.2 solar, respectively (Sakelliou et al 2002). This 
Mg/0 ratio derived from the RGS is a factor of 1.5 larger 
than that at r > 0.5' derived from the EMOS spectra 
using the APEC code for modeling the Fe-L. Since most 
of the photons detected by the RGS come from the very 
center, where a larger Mg/0 ratio is also allowed by the 
EMOS data, the Mg/0 ratio may increase at the center. 
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Fig. 13. The Mg/0 ratio derived from the line ratio, 
where the Fe-L lines are modeled with the APEC code. 
The meaning of the symbols are the same as in Figure 7. 
The dotted diamonds correspond to the abundance ratios 
of the EMOS derived from the APEC model fit to the 
deprojected spectra. Errors represent the 68% confidence 
level. 



5.6. Mg vs. Si 

The derived ratios of the Ka lines of H-like Mg and Si are 
summarized in Table 4 and Figure 14, Here, only the result 



using the APEC code for the Fe-L modeling is shown. The 
temperature dependence of the line ratio is quite similar 
to that of the Ka lines of O and Si, and the Mg/Si ratio 
was derived in the same way as in Sec. 5.1. Within 2', 
the contribution of the 1 keV component is taken into 
account. The derived Mg/Si ratio is ~ 0.6 solar, which 
agrees well with the ratio derived from the spectral fit 
using the APEC model in Sec 4. 

5.7. Fe vs. Si 

In contrast to the temperature dependences of the line ra- 
tios in the previous subsections, the ratio of the Fe-K line 
at 6.7 keV (which includes He and Li like ions) to the Ka 
line of the H-like Si line strongly depends on temperature. 
For example, at 2 keV, a 10 % increment of tempera- 
ture increases the Fe/Si ratio by 36 %. As a result, small 
uncertainties in the temperature structure give large un- 
certainties in the abundance ratios. Therefore, the Si/Fe 
ratio cannot obtained from the line ratio as in previous 
subsections. Since this ratio is quite important, we have 
calculated the Si/Fe line ratio of the projected data using 
the temperature structure obtained in Paper I and com- 
pared it with the observed profile (Figure |l^. We used 
the relation of kT = 1.68(1 + {R/1.6'f)°-^^keY, since 
the hotter component dominates these lines even within 
2'. Assuming the Fe/Si ratio is approximated by a -I- bR, 
where a and b are free parameters, we fitted the profile 
of the line ratio. The radial profile of the line ratio is well 
fitted with the model (x^=7.67 for 9 degrees of freedom) 
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Fig. 14. (Top panel) The ratio of strength of Ka hnes of 
H-hke Mg and Si plotted against the best fit temperature 
fitted with the single MEKAL model. The Fe-L is mod- 
eled with the APEC. Dashed lines correspond to constant 
abundance ratios. The meanings of the symbols are the 
same as in Figure 7. (Bottom panel) The radial profile of 
the Mg/Si ratio derived from the line ratio. The dotted 
diamonds are the Mg/Si ratios derived from the spectral 
fitting using the APEC code. Errors represent the 68% 
confidence level. 

and the Fe/Si ratio is determined to be 0.9 solar within 
the whole field of view (Figure |l^, although the abun- 
dance ratio changes by 20% when the whole temperature 
is shifted by 5%. Considering that the value of 0.9 solar 
is quite close to the ratio derived from the spectral fitting 
of the deprojected data, where the Fe abundance is deter- 
mined by the Fe-L emission, the Fe/Si ratio should be ~ 
0.9 solar and its radial gradient is less than 10-20%. 

6. Effect of resonance line scattering 

Some resonant lines may become optically thick in the 
dense cores of clusters. Shigeyama (1998) calculated the 
effect on M 87 and several X-ray luminous elliptical galax- 
ies, and found that it should be important at the center. 
Bohringer et al. (2001) suggested that the observed central 
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Fig. 15. (upper panel) The radial profile of projected line 
ratios of the EMOS (black) and EPN (gray) between Ka 
of H-like Si and He-like Fc. Errors correspond to 68% con- 
fidence level. The dashed line represents the best fit model 
using the temperature profile obtained in Paper I. (lower 
panel) 90 % confidence contour of the Fe/Si at i?=14' vs. 
that at R=l' obtained from the ratio of the Ka line of the 
He-fLi-like Fe to H-like Si of the projected data, assum- 
ing Fe/Si=a + bR and kT=1.68(l + (i?/1.6')^)°-"^ (keV). 
The cross is the correlation derived from the deprojected 
spectra, determined from the Fe-L emission. 



abundance drop may be due to the effect. In the case of 
an X-ray luminous elliptical, NGC 4636, Xu et al. (2002) 
discovered direct evidence of resonant scattering using a 
line ratio between optically thick and thin lines. 

The profile of optical depth for resonant lines of some 
prominent emission lines are calculated in Bohringer et al. 
(2001). Using the observed abundance profiles (Table 1), 
the temperature profile used in Sec. 5.7 and the density 
profile derived in Paper I, optical depths from the center 
for Ka lines of H-like Si and O are calculated to be 1 and 
0.4, respectively. Therefore, we must evaluate the effect 
of the scattering, although the central abundance drop is 



14 



Kyoko Matsushita et al.: XMM observation of M 87 



10 ^ F 




0.1 1 10 



r (arcmin) 

Fig. 16. Radial profile of brigtitness of H-like ions of O 
(bold black lines), Si (gray lines), and S (thin black lines). 
Crosses are observed profiles and solid and dotted lines 
corresponds to the simulated profiles before and after reso- 
nant line scattering, respectively. The bottom panel shows 
the ratios of the profiles of before and after the scattering. 

very small when fitted with the two temperature MEKAL 
model. 

Ignoring turbulent motion, we have calculated the ef- 
fect of resonant line scattering on the Ka lines of H-like 
ions, using Monte-Carlo simulations. 

The results are summarized in Figure The scat- 
tered profiles of lines are consistent with those ignoring 
scattering beyond 1'. Within this radius, the brightness of 
the lines decrease. The reduction rates at 0.5' are 10, 30, 
15% for O, Si, and S, respectively. The effect is important 
only within 0.5' on the Si line, where the uncertainty in 
temperature structure also gives an uncertainty in abun- 
dance determination. 

7. Discussion 

7.1. Summary of abundances 

The abundance profiles of O, Si, S, Ar, Ca and Fe are 
obtained from the deprojected spectra, based on the tem- 
perature study in Paper I. Although version 1.0 of the 
APEC code gives significantly different results (Paper I), 
version 1.1 of the APEC code gives consistent results with 
those of the MEKAL code for these elements. Within 2', 
considering the two temperature nature of the ICM, the 
abundance profiles from deprojected spectra become al- 
most constant. Within this radius, the Si, S, Ar, Ca and 
Fe abundances are ~ 1.5 solar, while the O abundance is 



a factor of 2 smaller. Using the APEC code, the Mg abun- 
dance is derived to be ~ 1 solar at the center, since the 
APEC code gives a better fit on the Fe-L/Mg-K structure 
at 1.4 keV. Beyond 2', Si, S, Ar, Ca and Fe show strong 
negative gradients, while the abundance gradient of O is 
smaller. The Si/Fe ratio is 1.1 solar, with no gradient in 
the field of view. 

The projected and deprojected profiles of abundance 
ratios are also obtained through line ratios. Since some 
line ratios among a-elements are nearly constant within 
the temperature range of the ICM around M 87, the abun- 
dance ratios are obtained with better accuracy. From the 
temperature dependence of the line ratio, we can conclude 
that the low 0/Si ratio cannot be explained by any tem- 
perature model, and the ratios of 0/Si are 0.4-0.5 solar 
at the center and 0.6-0.7 solar at the outer region. The 
Mg/0 ratio is determined to be 1.2-1.3 solar. The Ar and 
S abundances thus obtained are systematically larger than 
the spectral fit results due to small disagreements in the 
continuum fits around these lines. The S/Si is 1.1 solar at 
the center and 0.7 solar at outer regions. The Si/Fe ratio 
derived from the projected profile of the Ka line ratio is 
consistent with that derived from the deprojected spec- 
tra through the spectral fit, where the Fe abundance is 
determined by the Fe-L structure. 

The observed abundance pattern is similar to those 
obtained for A 496 (Tamura et al. 2001), and NGC 4636 
(Xu et al. 2002). Therefore, this abundance pattern may 
be uniform around the central galaxies of groups or clus- 
ters. 

7.2. Comparison of the abundance pattern with 
Galactic stars 

Figures |l^ and |l^ show the observed abundance ratios of 
O, Mg, Si, Ca and Fe compared to those of the Galactic 
disk stars (Edvardsson et al. 1993). The average abun- 
dance ratios of low metallicity (i.e. [Fe/II]<-0.8) stars of 
the Galaxy by Clementini et al. (1999), which should re- 
fiect the average abundance pattern of SN II products of 
our Galaxy are also plotted. Other papers on Galactic low 
metallicity stars give similar abundance ratios. [0/Fe] ra- 
tios for stars with [Fe/H]< —1 derived by Nissen et al. 
(1994) and Peimbert (1992) are 0.48±0.16, and 0.5, re- 
spectively. Gratton & Sneden (1991) gives [Si/Fe] of 0.3 
for these stars. 

The observed abundance pattern of M 87 is located at 
a simple extension of that of Galactic stars, although the 
observed Fe/0 range of M 87 is systematically larger. 

The Mg/0 ratio of the ICM is plotted assuming the 
ratio is constant within the observed region. We adopted 
the ratio derived from the line ratio, using the Fe-L mod- 
eling with the APEC code, since APEC gives better fits 
around the Mg-K/Fe-L region at 1.3-1.5 keV. [Mg/0] of 
the Galactic stars is almost constant at the same value 
as the ICM around M 87. The Galactic [Ca/Si] tends to 
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Fig. 17. The Fe/0 ratio is plotted against the Si/0 ratio. 
Those of disk stars (Edvardsson et al. 1993; closed circles) 
and average values of metal poor stars (Clementini et al. 
1999; open square) in the Galaxy are also plotted. The 
black solid line and dashed line represent the relation of 
the abundance pattern with (Fe/Si)sNia =1-1 solar (best 
fit relation) and (Fe/Si)sNia=2.6 solar (W7 ratio; Nomoto 
et al. 1984), respectively. The gray solid line and dashed 
line represent the relation of (Fe/Si)sNia=l-2 and 1.7 solar, 
respectively, adopting the abundance pattern of the SN II 
by Iwamoto et al (1999; asterisk) using the nucleosynthesis 
model of Nomoto et al. (1997) and assuming a Salpeter 
IMF. 
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a slightly smaller value than the that of M 87, but the 
difference is only 0.1 dex. 

We note that Xu et al. (2002) also discovered simi- 
lar values of the Fe/0 and Mg/0 ratio for the ISM in 
an eUiptical galaxy, NGC 4636. The Si/Fe ratio of this 
galaxy observed by ASCA is also determined to be 1 so- 
lar (Matsushita et al. 1997; 2000). Thus, the 0/Mg/Si/Fe 
pattern of NGC 4636 is similar to that of M 87. Therefore, 
the abundance pattern of the ICM of M 87 is not pecu- 
liar, but consistent, not only with an elliptical galaxy NGC 
4636, but also with that of our Galaxy. 

7.3. Nucleosynthesis from SN la 

The observed radial change of the 0/Fe ratio indicates 
that the contribution from SN la increases toward center. 

Although the abundance pattern of ejecta of SN II may 
differ between early-type and late-type galaxies, and that 
of SN la also may not be a constant (Umeda et al. 1999; 
Finoguenov et al. 2002), for a first attempt we have as- 
sumed (Fe/Si)sNia, (Si/0)sNii, and (Fe/0)sNii to be con- 
stants. Here, (Fe/Si)sNia is the Fe/Si ratio of ejecta of SN 
la, and (Si/0)sNii and (Fe/0)sNii are the Si/0 ratio and 
the Fe/0 ratio of the ejecta of SN II, respectively. Since 



[Fe/0] 

Fig. 18. [Mg/0], and [Ca/Si] are plotted against [Fe/0]. 
Those of disk stars (Edvardsson et al. 1993; closed circles) 
in the Galaxy are also plotted. 



O is not synthesized by SN la, the Fe/0 ratio is then 
expressed by. 



Fe _ /Fe\ 

O VSi/sNIa 



SNII 



Fc\ 

O / SNII 



(1) 



Thus, (Fe/Si)sNia is the inclination in the Si/0 vs. Fe/0 
plot (Figure |l^). Even within the M 87 data, the gradient 
indicates that (Fe/Si)sNia at the center is smaller than ~ 
1.3 solar. When we adopt the abundance pattern of the 
Galactic metal poor stars by Clementini et al. (1999) as 
that of SN II, (Fe/Si)sNia around M 87 is determined to 
be ~ 1.1 solar (Table ^ ). In addition, the M 87 data are 
located at the extension of those of Galactic stars, and 
the whole trend is consistent with (Fe/Si)sNia to be ~ 
1.1 solar. This means that also for metal rich stars in the 
Galaxy, the (Fe/Si)sNia is determined to be ~ 1.3 solar. 
These values are a factor of 2 smaller than the standard 
W7 model (Nomoto et al. 1984; Table ||), and closer to the 
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Fig. 19. The S/0 and Ar/0 ratio are plotted against the 
Si/0 ratio observed by the EMOS (black) and the EPN 
(gray). Dashed lines represents the average relation with 
90% significance. 



WDDl ratio (Iwamoto et al. 1999), which considers slow 
deflagration. 

In the same way, we have fitted the S and Ar to Si 
ratio (Figure |l9|) and derived the S/Si and Ar/Si ratio 
of SN la ((S/Si)sNia and (Ar/Si)sNia) to be ^ 1.5 solar 
and 1.3 solar, respectively (Table ^). In contrast to the 
Fe/Si ratio, W7 and WDDs give nearly the same values 
for the ratio between the intermediate elements, S/Si and 
Ar/Si, and the observed ratios are slightly larger than the 
theoretical values. 

7 A. Diversity of abundance pattern in SN la? 

The light curves of observed SN la are not identical but 
display a considerable variation (e.g. Hamuy et al. 1996). 
In SN la, the mass of synthesized Ni^^ determines the 
luminosity of each SN. Since the mass of the progenitor 



M 87 SNI 


W7 


WDDl 


WDD2 


WDDS 


1.57±0.32 


1.05 


1.14 


1.14 


1.14 


1.25±0.45 


0.76 


0.89 


0.91 


0.93 


1.10±0.25 


2.61 


1.34 


2.07 


2.99 



Table 5. The abundance pattern of SN la obtained from 
the ICM of M 87 and those from the nucleo synthesis 
model of W7 (Nomoto et al. 1984) and WDDs (Iwamoto 
et al. 1999). 



S/Si 
Ar/Si 
Fe/Si 



should be constant at 1.4 Mq, the ratio of mass of inter- 
mediate group elements from Si to Ca, to the mass of Fe 
and Ni, should depend on the luminosity of SN la. The ob- 
served luminosity of SN la correlates with the type of the 
host galaxy, and is suggested to be related to the age of 
the system; SNe la in old stellar system may have smaller 
luminosities (Iwanov et al. 2000), and hence are suggested 
to yield a smaller Fe/Si ratio (e.g. Umeda et al. 1999). 

As discussed in Finoguenov et al. (2002), the smaller 
Fe/Si ratio observed for the ICM around M 87 may re- 
flect the fact that M 87 is an old stellar system. When we 
compare the SN la abundance pattern for the central and 
outer regions of M 87, the outer zone SN la abundance 
pattern depends more heavily on the adopted abundance 
pattern of SN II. As in Finoguenov et al. (2002), we use the 
SN II pattern derived in Iwamoto et al. (1999) for which 
a Salpeter IMF was adopted. As indicated in Figure 17, 
the (Fe/Si)sNia of the outer region is determined to be 
^ 1.7 solar, while the central region data yield a value 
of ~ 1.2 solar (Figure |l^. In contrast, when we adopt 
the abundance pattern of metal poor Galactic stars as a 
SN II pattern, the difference of (Fe/Si)sNia of the central 
and outer regions becomes smaller. The average products 
of SN II may also differ between the ICM and late-type 
galaxies, since the IMF of stars may differ between early- 
type and late-type galaxies. Therefore, (Fe/Si)sNia of the 
outer region still has some uncertainty, while the central 
value is determined to be ~ 1.2 solar, which is not affected 
by the adopted SN II pattern very much. 

However, the metal poor stars, i.e. those with lower 
Fe/0, tend to be located around larger (Fe/Si)gNia values 
(Figure 17). This result suggests that the SN la were dom- 
inated by those with larger (Fe/Si)gNia when the Galaxy 
was a young stellar system. Furthermore, based on the 
ASCA survey of the Virgo cluster, Shibata et al. (2002) 
discovered that the Fe/Si ratio becomes larger at i? > 30', 
where the SN la products are accumulated from high z. 
Although the O abundance cannot be measured, this re- 
sult may reflect the fact that the Fe/Si ratio in SN la 
depends on the age of the system. 

The ICM in the central and outer regions should have 
a different origin of metals. As discussed in Paper I, con- 
sidering the gas mass and stellar mass loss rate, most of 
the Si and Fe at the center come from present SN la in 
M 87 in the last few Gyr. In contrast, in the outer re- 
gions, metals from SN la in M 87 are accumulated over a 
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longer time scale, and in addition, other galaxies should 
contribute to the metals. Since the SN la rate is expected 
to be much larger in the past (Renzini et al. 1993), the SN 
la ejecta which occurred in much younger systems should 
be important in the outer region. 

The indication that the abundance pattern of SN la in 
the outer region of M 87 features a larger (Fe/Si)sNia ratio 
compared to the inner region as found above, can nicely 
interpreted in line with of these previous findings. The 
larger (Fe/Si)sNia ratio originates from a younger stellar 
population, so SN la products that have been produced 
earlier in the history of M 87 and the Virgo cluster are 
more widely distributed in the ICM. 

7.5. Nucleo synthesis from SN II 

We can also constrain the abundance pattern of SN II in 
the ICM. 

Since Mg and O are not synthesized by SN la, the 
Mg/0 ratio reflects those of the products of SN II. The 
observed Mg/0 ratio of ^--^1.25 solar at r > 0.5' derived 
from the line ratio agrees well with that of the Galactic 
stars, although a higher Mg/0 ratio is also allowed within 
r < 0.5' due to the uncertainties in the temperature com- 
ponent below 1 keV. This value is also quite similar to the 
Mg/0 ratio of 1.3±0.2 solar of the ISM of an X-ray lu- 
minous elliptical galaxy, derived from the RGS data (Xu 
et al. 2002). The central Mg/0 reflects the stellar abun- 
dance ratio of M 87, while that of outer region is the ratio 
of of metals in the ICM which is an accumulation of old 
SN II products ejected from galaxies. Therefore, at least 
for the Mg/0 ratio, there is no observable difference in 
nucleosynthesis products of SN II between the elliptical 
galaxies, M 87 and NGC 4636, the ICM and the Galaxy. 

The observed S/Si ratio decreases when the 0/Si ratio 
increases (Figure 8, |9|). Considering that the observed 
Fe/Si ratio is constant, this indicates that S/Si ratio is 
anti-correlated to the 0/Fe ratio, that is the contribution 
from SN II/SN la. From equation (1), using S instead of 
Fe, the S/0 ratio in ejecta of SN II ((S/0)snii) is obtained 
from the observed relation in Figure y. The Si/0 ratio 
of metal poor stars by Clementini et al. (1999) is ~ 0.7 
solar, while the nucleo synthesis model by Iwamoto et al. 
(1999) gives the value of 0.9 solar. Therefore assuming 
that (Si/0)sNii is less than 1 solar, the extension of the 
observed values in Figure 19 gives (S/0)snii less than 0.5 
solar. In the same way, (Ar/0)sNii is derived to be less 
than 1 solar. 

The observed S abundances of Galactic stars are gen- 
erally consistent with the Si abundances (e.g. Chen et al. 
2002, Takeda-Hidai et al. 2002), although S lines are very 
week. In contrast, the observation of type II planetary 
nebulae yielded a S/0 and Ar/0 ratio to be 0.4±0.2 solar 
and l.libO.5 solar (Kwitter & Henry 2001). The nucleosyn- 
thesis model of SN II by Nomoto et al. (1997) assuming 
Salpeter's IMF gives the S/0 ratio of 0.6 solar, while an- 
other model by Woosley & Weaver (1995) gives a larger 



ratio of 1 solar. Thus, there may be some uncertainty in 
the S/0 ratio synthesized by SN II even in the Galaxy, 
and it is also not clear that the Galaxy and the ICM have 
the same ratio. The fact that the S/Si ratio of the ICM 
observed with ASCA generally decreases toward outer ra- 
dius (Finoguenov et al. 2000) supports the low S/Si ratio 
in SN II, although ASCA cannot constrain O abundance. 
Therefore, it is important to study outer regions of clusters 
with the XMM-Newton where SN II may become domi- 
nant (e.g. Fukazawa et al. 2000, Finoguenov et al. 2000). 

In summary, the SN II abundance pattern observed 
from the ICM mostly agrees well with that of our Galaxy. 
Since the abundance pattern of SN II depends on the IMF 
of stars, this consistency is able to provide strong con- 
straints on it. 

7. 6. Comparison with the stellar metallicity profile 

Figure |2^ shows the observed Mg profile, with the stel- 
lar metallicity profile from the Mg2 index (Kobayashi & 
Arimoto 1999). In addition to the Mg abundance profile 
derived from the APEC model fit, that from the Mg/Si 
line ratio and deprojected Si abundance profile is also plot- 
ted. The Mg abundance profile of the ICM has a slightly 
smaller normalization by 20~30% than the stellar metal- 
licity profile at the same radius. 

Although the Mg abundance may have some uncer- 
tainties due to the Fe-L atomic data, the APEC code 
can well fit the Fe-L/Mg-K structure. The uncertainty of 
the Mg abundance due to the uncertainty of the tempera- 
ture structure should be also small, at least at i?=0.5-2', 
since the three temperature model does not change the 
result, although the abundances from the two- and single- 
temperature results differ by a factor of 1.5-2. Considering 
these uncertainties and the difference in observational 
techniques, we can conclude that the Mg abundance of 
the ICM is consistent with the stellar metallicity profile 
at the same radius. Since we are comparing abundances in 
two distinct media, stars and ISM, which could have very 
different histories, the abundance results do not have to 
agree in general. But this agreement is consistent with the 
picture where the central gas in M 87 comes from stellar 
mass loss as discussed below. 

7.7. The SN la rate of M 87 

The metals in the ICM are a sum of metals ejected from 
the galaxy and those contained previously within the 
ICM. The observed Fe abundance becomes nearly con- 
stant within 2', which indicates that the gas is domi- 
nated by gas ejected from the galaxy. The observed Mg 
and Fe abundances within 0.35-2' are 1.1 ± 0.2 solar 
and 1.54±0.1 solar, respectively. Both the nucleosynthesis 
model by Nomoto et al. (1997) and observation of Galactic 
stars (e.g. Clementini et al. 1999; Nissen et al. 1994) indi- 
cate that (Fe/Mg)sNii is 0.3-0.5 solar. Therefore, within 
R < 2' , the Fe abundance synthesized by SN II should 
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Fig. 20. Mg abundance profile of the EMOS data derived 
from the spectral fitting with the APEC code (diamonds) 

and from line ratios (crosses). The dashed line represents 
the stellar metallicity derived from Mg2 index (Kobayashi 
& Arimoto 2000) 



theoretical stellar evolutionally model for a stellar pop- 
ulation, the stellar mass loss rate is approximated by 
1.5 X 10~^^ LBti^'^MQ/yr, where ti^ is the age in unit 
of 15 Gyr, and Lb is the B-band luminosity (Ciotti et al. 
1991). An infrared observation, which can directly trace 
mass loss rate from asymptotic giant branch stars, gives 
a consistent value within a factor of 2 with the theo- 
retical value (Athey et al. 2002). The SN la rate of el- 
liptical galaxies is estimated to be 0.13±0.05 /lys^ SN 
Ia/lOOyr/lOi°L0 (Cappellaro et al. 1997). Here, hrs is 
a Hubble constant in unit of 75 km s~^ Mpc~^. Using the 
Hubble constant of 72±8 km Mpc"'^ (Freedman et 
al. 2001), the present contribution from SN la to the Fe 
abundance becomes 1.4-3.6 solar. 

Considering the large error due to a statistical error 

and uncertainties in bias corrections of SN la from optical 
observation, this value is consistent with the Fe abundance 
synthesized by SN la within 2', 0.9-1.4 solar, which is a 
sum of ejecta of present SN la and that was trapped in 
stars and comes from stellar mass loss. The contribution 
from stars to Fe synthesized by SN la may be small, since 
stars in giant elliptical galaxies are thought to have a SN 
II like abundance pattern (e.g. Worthey et al. 1993). 



be 0.3-0.6 solar. Subtracting it, the Fe abundance synthe- 
sized by SN la is derived to be ~ 0.9-1.4 solar. 

The metallicity of gas ejected from the galaxy are a 
sum of stellar metallicity and the contribution from SN 

la, and can be expressed as. 



z = 



a* 



■ y* + (— >ysN 



a* 



(2) 



(Loewenstein and Mathews 1991; Ciotti et al. 1991), where 
is the mass fraction of the ith element, a* and asN are 
the mass loss rates of stars and SNe respectively, and yX 
and i/gj^ are their yields, is the SN la rate and Mgj^ 
is the mass of the zth element synthesized by one SN la. 
Converting to abundance in unit of the solar ratio, Fe 
abundance synthesized by SN la of gas ejected from an 
elliptical galaxy becomes 



^* SNIa 



4Fe 

^* SNIa 



xlO-i3/yr/I,B 



V solar } 

\ ( <N \ 



5xlO-iiMo/yr/LE 



(3) 



Here, ^^''gNia is stellar Fe abundance synthesized by 
SN la and Mg^ is the mass of Fe synthesized by one SN 

is the Fe metallicity of gas with the solar abun- 



la. 



^Fo 
''solar 



dance. Considering the fact that recent observations of 
M 87 and A 496 (Tamura et al. 2001) indicate a smaller 

Fe/Si ratio synthesized by SN la than usually assumed 
previously, Mpe of one SN la should be ~ O.4M0. From a 



7.8. Implication of S abundance 

In Sec. 7. 5 we indicated that the S/Si ratio in SN II is much 
lower than 1 solar. If this fact is valid universally in the 
ICM, then the S/Si ratio may be a better indicator for the 
contribution from the SN la/SN II than the Fe/Si ratio. 
Of course, the best indicators of SN II are O, Ne and Mg, 
which are not synthesized in SN la. However, emission 
lines of Ne and Mg are often hidden in the strong Fe- 
L lines, and the O abundance can be measured only in 
relatively cool clusters. 

As discussed in Sec. 7.4 the Fe/Si ratio of SN la may 
be expected to show a variation, since the luminosity of 
the SN la should reflect the amount of Ni'"'^ synthesized 
by SN la. The pure SN II ejecta may provide the smaller 
Fe/Si ratio. Then SN la in young systems occur, which 
may have a larger Fe/Si ratio and so the Fe/Si ratio in- 
creases. Finally, SN la in old systems decrease the Fe/Si 
ratio again. In contrast, at least in Iwamoto et al. (1999), 
the S/Si ratio in the various SN la is almost constant. In 
addition, as shown in Sec. 5.2, the S/Si line ratio does 
not depend on the plasma temperature, while that of the 
Fe/Si is a very strong function. This means that, the de- 
rived S/Si does not depend on the temperature structure 
of the ICM. This is a very important advantage for the 
cool cores of the cluster center, where many temperature 
components exists. 

The problem is that both Si and S are synthesized in 
the same nuclco process, i. e. in the same region in the 
SN. Therefore, it is surprising that the S/Si ratio shows 
a radial variation. Unfortunately, the S abundance mea- 
sured of the Galactic stars has very large uncertainties, 
and the S/Si ratio of the SN II may also depend on the 
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IMF of stars. Therefore, it is quite important to cahbrate 
the S/Si vs. 0/Fe pattern in luminous clusters. 

8. Summary and Conclusion 

Based on the temperature structure studied in Paper I, 
abundance profiles of 7 elements of the ICM around M 87 
are obtained using the deprojected spectra. We have dis- 
cussed the use of one- and two-temperature models in the 
fit and the application of the MEKAL and APEC codes. 
Two temperature models are only important in the re- 
gion, R< 2 arcmin. The previous problems with an earlier 
version of the APEC code are now resolved in the present 
version and this code provides now better fitting solutions 
in general. In addition, using line ratio profiles of the pro- 
jected and deprojected data, the abundance ratio profiles 
are obtained with smaller uncertainties, considering the 
temperature dependence of the line ratio. The main re- 
sults obtained are as follows, 

— Abundance profiles of 7 elements are obtained with 
high accuracy (e.g. ,5Fe~ (5Si~5%, (5O~10%). 

— The Si and Fc abundance profiles have steep gradients 
at > 2' with a Fe/Si ratio of 0.9 ±0.1 solar, and become 
constant within this radius at 1.5-1.7 solar. 

— The S/Si ratio is about 1 solar at r < 4' and decreases 
to ~ 0.7 solar at r > 10'. 

— The Ar/Si and Ca/Si are about 0.8 solar and 1.5 solar, 
respectively, although errors are relatively large. 

— The O and Mg abundances are smaller than the abun- 
dances of Fe and the intermediate elements. The 0/Si 
ratio is less than half solar at the center and increases 
with radius. The Mg/0 ratio is 1.25 solar and consis- 
tent with no radial gradients at least at r > 0.5'. 

— The observed abundance pattern among O, Mg, Si, Ca 
and Fe of the ICM is located at an extension of that of 
Galactic stars. Therefore, the abundance pattern of the 
ICM is not peculiar and it should strongly constrain 
the products of SN la and SN II of early-type and 
late-type galaxies. 

— The Si/O/Fe diagnostics shows large Si contributions 
by SN la, which may reflect the observational finding 
that SN la in old stellar systems are fainter. 

— The observed Mg abundance of the ICM is consistent 
with the stellar metallicity profile from the Mg2 index 
at the same radius. This result is consistent with stellar 
mass loss as the source of the gas in the very central 
region of M 87. 

— The central abundance can be explained with the ob- 
served SN la rate and SN la model yields. 

— The different radial profiles between Si and S suggest 
that the S/Si ratio synthesized from SN II is much 
smaller than the solar ratio. Then, the S/Si ratio may 
be a better indicator of the relative contribution from 
SN la and SN II, when O and Mg abundances cannot 
be observed. 

We are actually very fortunate to have M 87 and the 
Virgo cluster as the closest galaxy cluster in our neighbor- 



hood, since it offers an ICM in just this low temperature 
range where the spectra are richest in spectral lines for 
this type of abundance diagnostics. Therefore it would be 
worthwhile to use this unique case for even deeper obser- 
vations and more detailed investigations by X-ray spec- 
troscopy in the future. 
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